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Abstract 

We introduce a semiclassical quantization method which is based on a stro- 
boscopic description of the classical and the quantum flows. We show that 
this approach emerges naturally when one is interested in extracting the en- 
ergy spectrum within a prescribed and finite energy interval. The resulting 
semiclassical expression involves a finite number of periodic orbits whose en- 
ergies are in the considered interval. Higher order corrections which reflect 
the sharp restriction of the spectrum to an interval are explicitely given. The 
relation to Fourier methods for extracting semiclassical spectra, such as har- 
monic inversion, is worked out. The constraints due to the finite dimension 
of the Hilbert space and the unitarity of the restricted quantum evolution 
operator are important ingredients in this context. 
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I. INTRODUCTION 



The pioneering semiclassical quantization schemes introduced by Gutzwiller Jl|J^] (general 
chaotic hamiltonian flows) and Balian and Bloch (billiards) || are afflicted by a severe and 
intrisic problem: the semiclassical trace formulae and the spectral determinants derived from 
them involve an unlimited number of periodic orbits and do not converge on the real energy 
axis. Several authors proposed various solutions to this problem: 

The cycle expansion j§-§| exploits the presence of correlations between long and short 
orbits to argue that a suitable rearrangement of the spectral determinant converges on the 
real axis. The improved convergence has been demonstrated for certain open systems like 
scattering off three disks |||7| and for a map, the quantum bakers map ||. Applications to 
bounded systems are typically troubled by the lack of a convenient symbolic coding and the 
presence of marginally stable periodic orbits that spoil the ideal exponential convergence 

MM- 



Keating and Berry fllT|-|i"3|j made use of the analytic properties of the exact spectral 
determinant, and enforced it on the semiclassical expression, thus deriving a quantization 
scheme which is similar to the Riemann-Siegel formula for the Riemann ( function on the 
critical line. Its main term involves a finite number of periodic orbits, and the further 
corrections are given explicitly and they are small in the semiclassical limit. 

There are, however, alternative semiclassical quantization schemes for chaotic systems: 
One considers the quantum analogue of the classical evolution operator on a Poincare surface 



of section of a finite volume [0-|T6[. The area preserving Poincare map when quantized 



yields a semiclassically unitary evolution operator which acts on a Hilbert space of dimension 
L = , where V is the phase space volume of the 2f dimensional section, and [•] stands for 
the integer part. Denoting by S(E, h) the Lx L (semiclassically) unitary evolution operator, 
one expresses the quantization condition in terms of the secular function: 

Z sd (E;h) = det(I-S(E,h)) , (1.1) 

and the spectrum {E n } satisfies the quantization condition 

Z scl (E n ,h) = . (1.2) 

The secular function is nothing but the characteristic polynomial of S, det(I — zS(E, h)) = 
Y J f=Q a i{E. l h)z l computed at z = 1. The special symmetry of the Riemann-Siegel expression, 
which is a consequence of the functional equation in the Riemann-^ function case, follows 
here from the unitarity of 5*. It implies the inversive identities 

a x = e i0 a^ , (1.3) 

where G stands for the phase of det (—S(E,h)). When the semiclassical approximation is 
used to compute Z sc i(E; h), the coefficients ai are expressed in terms of periodic orbits of 
the Poincare map with periods up to I. Thus, the computation of the semiclassical secular 
equation requires only a finite number of periodic orbits. 

This proceedure was used e.g. by Bogomolny who derived a semiclassical approx- 



imation for the restriction of the Green function to the Poincare section, and by Doron 
and Smilansky [Rj| who used the scattering approach, and identified the scattering matrix 
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as the unitary operator in equation ( [L.l]) . Further work along this line was carried out by 
Gutzwiller [[Hj, Prosen [[L8[] and others (for a review see e.g., jlBj). 

The reduction of the classical dynamics to a discrete, area preserving mapping, can also 
be achieved by observing the flow at fixed time intervals. It is our aim here to suggest 
a new quantization scheme which is based on the quantum mechanical analogue of the 
stroboscopic classical map. The resulting quantization scheme, which is developed in section 
(y), addresses the spectrum in a finite spectral interval of a size which is dictated by the 
dynamical system and by Planck's constant. It is similar in its formal structure to the 
methods which were briefly summerized above. However, it is based on the classical periodic 
orbits of a different map, and it offers both the leading and the next to leading order 
contributions in the semiclassical expression for the secular equation. The semiclassical 
version of the stroboscopic quantization is given in (||). Among the other consequences 
of our derivation we note a relation between the time steps in the evolution operator and 
the mean density of states, a consistent truncation of the periodic orbit series as well as 
conditions for self-inversiveness of the characteristic polynominal. 

The idea of extracting eigenvalues within a certain energy band has many connections 
to frequently employed numerical schemes that also obtain eigenvalues within a certain fre- 
quency band. For instance, from the overlap of a wave packet with its propagated image 



over some time interval one can extract eigenvalues in a band near the initial energy [19 



The semiclassical procedure described here achieves the same without a weighting by the 
projection of initial wave packet onto eigenstates and moreover suggests a semiclassical ap- 
proximation. A numerically accurate and efficient way of implementing this Fourier analysis 
has been developped by Neuheuser and Wall pO| , |2T|j and Mandelshtam and Taylor [p2]-|2""l| 
in the form of harmonic analysis. The way harmonic analysis beats the resolution limit set 
by the Nyquist theorem is by implementing additional information on the Fourier signal. In 
the case of the traces it uses explicitely that they are sums of exponentials with a discrete 
set of frequencies. We will relate the basic ideas of harmonic analysis to the semiclassical 
traces and discuss the relation to the cycle expansion. The relation of the present formalism 
and harmonic inversion will be discussed in section (|TV|) . A few concluding remarks and a 
summary will be the contents of section (|V|) . 



II. STROBOSCOPIC QUANTIZATION 

We would like to compute a finite sequence of energy eigenvalues in the spectrum of an 
autonomous (time independent) quantum Hamiltonian H. This sequence is located in an 
energy interval S, 

e = [Eo- A/2 , E + A/2] . (2.1) 

The midpoint of the interval Eq is arbitrary, and its length A is ~ (9(/i 7 ) where 1 > 7 > 0. 
This makes A small on the classical scale but large on the quantum scale since for a system 
with / degrees of freedom the mean spectral density d(Eo) ~ 0(h*). Hence the number of 
eigenenergies in the interval 

N = d(E )A, (2.2) 
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is large. Moreover, A is also sufficiently small so that d(E) is effectively constant in the 
interval. 

Consider now the N dimensional Hilbert space spanned by the eigenvectors of H with 
eigenenergies in £. The quantum evolution operator in this subspace can be projected from 
the full evolution operator U(t) = e - lHt / h , it is given explicitly by 

U(t;S) = £ \^n)e- iE ^ h (^ n \ (2.3) 
E n e£ 

and 

txU{t;£) = e ~ lEnt/h = [ X{E) d(E)e' iEt / h dE , (2.4) 
where x{E) is the characteristic function of £, and 

oo 

d(E) = J2 5 ( E ~ E ") ( 2 - 5 ) 

v=l 

is the density of states for the full system. We will use small greek indices to run over 
all eigenstates of the Hamiltonian and small latin ones if they are restricted to the energy 
interval £. 

Choose an as yet arbitrary time interval of length r so that r A < IitTi. The N eigenvalues 
e -iE„r/h Q f foe unitary matrix U{r\ £) occupy an arc of the unit circle whose length is at 
most 27r, and therefore they preserve the order of the spectral sequence on £. The spectrum 
of H in the interval £ can be identified as the zeros of the secular equation: 

( £ (E)=det(l-e^U(T;£)) . (2.6) 

In other words, whenever the eigenphases of e% Er U{r\ £) take the value 1, a spectral point 
of H is encountered. Because of the exponentiation and the restriction to a finite set of 
eigenenergies in the interval £, the density of states derived from this secular equation is a 
train of delta functions at the positions of the eigenenergies E n , periodically continued by a 
separation h/27rr, 

oo 

d(E; £) =r/h £ £ 5{ r{E - E n )/% + 2vrm) 

m=-oo n:{E n e£ } 

= ^+^ £ { e-.-/» tlU{mTl S) + c,.} . 

The second line follows from Poisson summation of the first line. As usual, the density 
of states can be split into a smooth part and an oscillatory part, where the latter should 
vanish when averaged over the energy interval H/27TT. In the semiclassical limit, the traces 
trW m (r; £) are expressed as sums over periodic orbits, which give the oscillatory part of the 
spectral density. The smooth parts of both sides of ( |2.7|) must be identical, which defines r 
uniquely as 
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27rhd(E ) T H 

T = ^r~ = iv' (2 - 7) 

where th is the Heisenberg time. 

This choice of the value of r completely determines the quantum map U(t; £) which 
provides us through (|2.7|) with the spectral density of H in £. This is a stroboscopic map, 
since r is fixed. The unitarity of U{r\ £) is an important asset, since one can write the secular 
equation in a form which exploits naturally the inversive symmetry of the coefficients of the 
characteristic polynomial of U[t\£\ Writing 



N 

Ce(E) = det (/ - e lEr / h U(r; £))='£ aJ nEr ' n , (2.8) 

n=0 



the inversive symmetry is expressed by the identity 

a n = e l0 a* N _ n , (2.9) 
where stands for the phase of det(— U(t; £)). We use 

det (-jEr/n u[r . £ ^ = e ^N + i2n(E-E )d(E ) = ^N(E,E ) ^ ^ 1Q) 

where we introduced the smooth spectral counting function in the energy interval £, 
N(E, E ) = {E- E )d(E ) - N/2. Thus, 

[AT/2] 

e i7r £ (E) = {«ne i7r N{E ' Eo) + a* n e-™ N ( E > E °)} + e N a N/2 e™ N ^ E ^ . (2.11) 

n=0 

As before, [x] stands for the largest integer smaller than x, and = 1 if is even and 
ejy = otherwise. 

The cycle expansion [p6|-p8| and Newton's identities relate the traces t n = tilA n {r;£) 
and the coefficients of the characteristic polynomial 

a n = [t n + ^2a k t n _ k \ . (2.12) 

n \ k=i ) 

The inversive symmetry halves the number of coefficients needed to calculate the spectral 
secular equation. This symmetry is a consequence of the unitarity of the evolution operator 
and as such it is the expression of a basic property of the quantum evolution - the conser- 
vation of probability (norm). The semiclassical theory provides approximate expressions to 
t n which are used to compute the a n using Newton's identities. The resulting semiclassical 
secular function is similar in many ways to the expressions derived by Berry and Keating 
TTHl3[, Bogomolny and Doron and Smilansky [p~5],p~6|| . It will be discussed in the sequel. 



III. THE SEMICLASSICAL APPROXIMATION 

It is convenient to introduce the semiclassical approximation by writing (|2.4| ) as 
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txU{t-£) = J x (E)d(E)e- iEt/h dE 
1 r°° 

-r / ds x(s)trU(t + s) 
Inn J-oo 



where 



sin SA 

X (E)e* Es dE = e iEos / h =^2L 

2h 



(3.1) 



(3.2) 



and where U(t) = e lHt / h is the evolution operator of the entire system. We can then use 
the semiclassical expression for tr U (t) for any time t, given by 



[tr U(t)} 



scl 



t I dE \ 1 ^ QiRp/H-iv P K/2 



(27T^)V2 



P ev(t) 



\dt J \det(I-M p 



1/2 



(3.3) 



where V(t) is the set of t-periodic orbits and 



Rrt 



p p (t)q p (t)-H(q p (t),p p (t))]dt . 



(3.4) 



The building blocks for the semiclassical theory are the traces trU n (r; £). By an exact 
quantum identity they coincide with trU(nr;£), so that we need to compute the traces 
for times which are integer multiples of r. Notice that the periodic orbits p G V(t) which 
contribute to (3.3) can have arbitrary energies E p (t). They are functions of the period t = nr 



which is fixed. Substituting ( |3~3] ) in ( |3.1[ ), 

1 r°° 

[trU(nr;£)] sd = — / ds X {s) [trC/(nr + s)} t 



2/1 II J-oo 
j roo 



2nh 



/oo 
ds x(s) 
-oo 



/ &E \ 1,/2 e iRp/h-iu p Tr/2 



nr + s . 
(W 71 P evh +S) W \det(I-M p y 12 



(3.5) 



The range of x(s) is determined by the time ^ ~ (9(/i 1_7 ) which is classically short. Hence 
we can approximate the integral over the periodic orbit sum by expanding the phase of the 
integrand about t. For this purpose we use 9R ^ = —E p {t). Denoting the preexponential 
amplitudes in ( p.5| ) by A p (t), we use A p (nr + s) « A p (nr) which is consistent with the 
semiclassical approximation. The contribution of each periodic orbit will be 



2nh 



(E -E p (t))s-^ S 2 



2h 



1 

-i \ 2 



A p (nr) e^(-) {[C(x+(nr)) + iS(x+(nr))] - [C(xj(nr)) + iS^M)] } , (3.6) 



where C and 5 are the cos and sin Fresnel integrals with the argument 

± Sp(nr) - ± f 
(nr) = — ^ 



(2ft 



dt 



(3.7) 
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In the above expression, the difference of the energy of the periodic orbit, E p (nr), from 
the interval of interest 8, is measured in units of the energy scale (2h 



2. (we assumed 

throughout the derivation that, as is typically the case, is negative.) Using dimensional 

analysis one can bound this scale from above by O (fi^ 1 ^, which is small on the classical 
scale. Using the asymptotic expression for the Fresnel integral for large argument, 

C(x) + iS(x) « (-Y sign(x) - — e^ x2 + O ( ^V) • (3-8) 
\2J xti \ \xr I 



The difference between the Fresnel integrals in ( |3.6| ) is dominated in the semiclassical limit 
by the contribution of the leading term in ( |3.8j ) which, together with the factor (— i/2) 1 ^ 2 (see 
( p.6[) ) is the characteristic function of the spectral interval 8. Thus, the leading semiclassical 
contribution makes use of periodic orbits whose energy E p (t) is in 8 and does not take into 
account the "diffractive" effects due to the sharp restriction of the spectrum to the interval 
8. The Fresnel functions with finite arguments include the appropriate corrections. 

The semiclassical approximation for the [t n ] sd = [tr U(nr; 8)] sd is obtained by summing 
the amplitudes Q3.6|) over the set V(nr) of periodic orbits. When [t n ] sd is substituted in 
( |2.11| ) one obtains the semiclassical spectral secular equation. It can be written as a sum of 
two parts. The first is obtained from the leading semiclassical expression for t n . It uses the 
periodic classical orbits of period nr with energies in the interval of interest 8. The diffractive 
corrections go beyond this limit, and apart from modifying the contributions of "allowed" 
periodic orbits, it introduce the effects of periodic orbits whose energies are outside the 
strict energy interval. These corrections are analogous to the expressions derived by Berry 
and Keating [|Tl~1 -|T3|j which are missing in the derivations based on the Poincare section 



or the scattering approaches to quantization [ 13 ,|TB | . Here, one can attribute them to the 
sharp truncation of the energy domain and their main effect is the inclusion of evanescent, 
classically forbidden contributions. 

Finally, we would like to mention that the semiclassical approximation could be intro- 
duced in a different way. With 

P{£)= £ \n)(n\, (3.9) 

n:{E n e£} 

the projector onto the energy interval, we have the quantum identity 

trW n (r; 8) = trW(nr; 8) = tr {P{8)U{r)) n . (3.10) 
The semiclassical approximation we used in (|3~3|) and ( |3.5|) is 

trW»(r) ~ [tTU(nr)] scl . (3.11) 

The alternative approximation 

trW»(T) ~ tr ([P(8)U(r)] scl T . (3.12) 

would contain a product of Fresnel factors that reduces to ( |3.5| ) only to leading order. Such 
differences are not uncommon in attempts to go beyond the leading order in the semiclassical 
approximation. We favor the representation ( |3.11| ) since it makes use of as much exact 
quantum information as possible and leaves the semiclassical approximation only to the 
very end. 
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IV. OBTAINING EIGENVALUES 



Given the traces of U(nr) at equidistant time intervals r, the extraction of eigenvalues 
becomes a problem in Fourier inversion. If all traces tiU m in ( |2.7| ) and ([2.12| ) are known, the 



relation is exact. However, in practice only a finite number of traces can be calculated, and 
then standard Fourier inversion is limited in resolution by the Nyquist sampling theorem, 
which for traces up to N implies a resolution of the order of the mean spacing. Additional 
errors are introduced by the semiclassical approximation. In order to go beyond the limit of 
resolution set by the Nyquist theorem additional information has to be added to the Fourier 
inversion. For the case at hand, harmonic inversion |2D|-F25J is particularly appropriate. It 
assumes that the signal is composed of a finite number of frequencies and provides an efficient 
method for solving the ensuing approximation problem. Specifically, it uses the traces to 
set up a matrix V that has the same eigenvalues as U, so that the frequencies can be found 
from the traces by an eigenvalue determination, which numerically is more reliable than a 
standard search for zeroes. Note that the matrix V has the same eigenvalues as U, but it 
will usually not be unitary. In the description of harmonic analysis below we will emphasise 
the matrix structure and its origin and will omit numerical issues such as windowing and 
the like. 

Let iV be the fixed dimension of the matrix U and let, as before, 

N 

t m = tiU m = ]Te-* m ^, (4.1) 

n=l 

be the traces of the m-th power of the unitary operator. In the quantum case the phases 
<f) n = E n r/h are real and contain the eigenenergies. In the semiclassical approximation they 
may become complex, in which case the imaginary parts can be used as a measure of the 



semiclassical error. The form ( f4.1| ) of the traces implies that t n for n > N can be expressed 
as a linear combination of the traces t m with indices m < N, as in the case of autoregressive 
models. Harmonic inversion exploits these relations for the construction of a matrix V that 
has the same eigenvalues as U. To derive this matrix, define an iV-dimensional vector t m of 
traces starting with t m , i.e., 



(4.2) 



\tm+N-l 

the iV-dimensional vector 1= (1, • - • , 1) T and the auxiliary N x N matrix 

On,„ = e- im4>n ■ (4.3) 

Then 

ti = fil. (4.4) 
For higher m the vectors t m can be represented as 

t m = nV m ~ 1 I. (4.5) 



s 



with the diagonal matrix 



V = e 

' m,n ° 



-i<t>nX 



(4.6) 



The aim is to find a matrix V that is similar to V, so that the frequencies can be determined 
from the eigenvalues of V. To this end solve (|4.4|) for I and substitute in (|4.5|) : 



(4.7) 



Since also 



t m+k = ov m - 1 v k n~% = (pv^or 1 ) (nv k orHi 



k+l 



it is possible to construct the full matrix 



(4.8) 



(4.9) 



from the images of the vectors t n for different n. Specifically, for m — 2, the entries for the 
first iterate of V are 



k+l 



(4.10) 



This result shows that the matrix V\ can be constructed from the traces (assuming that the 
eigenvalues are not identical and the vectors of traces not linearly dependent). The matrix 
Vi is neither unitary nor symmetric but nevertheless has the eigenvalues exp(— i<p n ). 

It is possible to proceed one step further and to derive the characteristic polynominal. 
With the N x N matrices 



T. 



tM+l 



M 



t]M-l 



\tM+N-l tM+N-2 



tjM-N+1 \ 
thd-N+2 



(4.11) 



the relation ( |4. 10|) becomes 



T 



N+l 



ViT, 



N 



(4.12) 



so that formally V\ = Tjv(Ijv+i) 1 . Since — 1 columns in T/v+i and Tn coincide (up to a 
shift to the right), V\ has the form 



Vx 



( 





1 








1 








1 



\b N fejv-i bN-2 b N -s 
where the vector b = (pN, &at-i, • " ,bi) solves 



0\ 



bij 



(4.13) 
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( t2N-l \ 
t2N-2 

V t N -i J 



■TV 



( cln \ 

O-N-l 



\ «1 / 



(4.14) 



Exploiting this special form of the matrix V\ and the vector b the Fredholm determinant 
can be expanded 



N 



F(z)=det(l-*Vi) = 



(4.15) 



k=l 



Since the eigenvalues of U and V\ are the same, the coefficients of the characteristic poly- 
nomials (|2.8|) and (|4.15|) also have to be the same, up to an overall scale factor. 

It is important to note one significant difference, though: in the derivation of the poly- 
nominal for the quantum operator U, traces up to order N determined all coefficients and 
self inversivness allowed to bring this number down to N/2. In the harmonic inversion case 
the number of traces needed is 2N+1. The origin of this difference is the fact that in the case 
of the harmonic inversion nothing changes in the formalism if the traces trU n are replaced 
by trAU n with some operator A: the set up for the matrix V\ remains the same, but the 
coefficients of the vector I change. If A = 1, then traces with n > N can be expressed by 
traces of lower order using Cayley's theorem. Thus, harmonic inversion is more general in 
its basic structure, but also more wasteful in the number of traces required. In addition, it 
is difficult to see how properties such as self-inversiveness can be implemented directly (they 
can always be put in by hand in the characteristic polynominal) . 



V. CONCLUDING REMARKS 

The main features of stroboscopic quantization are the limitation to a finite interval in 
energy, the semiclassical expression with a finite number of periodic orbits and the finite 
characteristic polynominal. The origin of the finite representation is a physical one, the 
limitation to a finite interval in energy. Other modifications of the Gutzwiller trace formula, 
such as a Gaussian truncation or even a sharp cut-off, arrive at this restriction at the price 
of smearing the eigenvalues. Moreover, since with these smearings the evolution operator is 
not confined to a finite interval, the characteristic function is not a polynominal anymore 
and the conditions of self-inversiveness are difficult to implement. 

The weight of periodic orbits and in particular the Fresnel corrections bear a striking 
similarity to the error function truncations introduced by Berry and Keating on the level 
of the Fredholm product fTTHD! • Note, however, that the product of two Fresnel integrals 



is not a Fresnel integral of the sum of the arguments, so that the modification that would 
apply to a pseudo-orbit differs from the product of the weights of the original orbits. 

On the numerical side, the results on stroboscopic quantization suggest that for harmonic 
analysis the time period r over which the evolution operator is followed should be fixed to 
be th/N, so that an optimal characteristic polynominal that shows self-inversiveness results. 
Ideally, one would like to improve the method so as to use only the linearly independent 
N/2 traces, but it is not clear how that can be implemented other than by brute force. 
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The finiteness of the characteristic polynominal results from the quantum fact that after 
the projection onto the energy interval the Hilbert space is finite dimensional and that there 
are only a finite number of eigenvalues. In particular, this implies relations between traces 
of higher powers of U due to Cayley's theorem. Within the semiclassical approximation 
no such constraint and no such relation between traces of higher and lower powers are 
evident. In many approaches [|TT| p~6]| they are put in by hand and justified by appeal to 
the properties of the quantum propagator. However, it is legitimate to ask to which extend 
the semiclassical approximation reflects these properties by itself, and it is comforting to 
note that, as demonstrated for the Bakers map || and also for the three disk scattering 
system PJ7|J30[], where none of these unitarity arguments apply, the coefficients in the cycle 
expansion beyond the dimension of the system decay. 
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